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Abstract. Studies of distant galaxies have shown that ellipticals and large spirals (Schade et al. 1999; Lilly et al. 1998) 
were already in place 8 Gyr ago, leading to a very modest recent star formation (Brinchmann and Ellis, 2000) in intermediate 
mass galaxies (3-30 10'" Mq). This is challenged by a recent analysis (Heavens et al. 2004) of the fossil record of the stellar 
populations of ~ lO"" nearby galaxies, which shows that intermediate mass galaxies formed or assembled the bulk of their stars 
4 to 8 Gyr ago. Here we present direct observational evidence supporting this findings from a long term, multi-wavelength 
study of 195 z> 0.4 intermediate mass galaxies, mostly selected from the Canada France Redshift Survey (CFRS). We show 
that recent and efficient star formation is revealed at IR wavelengths since ~ 15% of intermediate mass galaxies at z> 0.4 are 
indeed luminous IR galaxies (LIRGs), a phenomenon far more common than in the local Universe. The star formation in LIRGs 
is sufficient in itself to produce 38% of the total stellar mass of intermediate mass galaxies and then to account for most of the 
reported stellar mass formation since z=l. 

Observations of distant galaxies have also the potential to resolve their star formation and mass assembly histories. The high 
occurrence of LIRGs is easily understood only if they correspond to episodic peaks of star formation, during which galaxies 
are reddened through short IREs (infrared episodes). We estimate that each galaxy should experience 4 to 5 x (tjre/O. IGyr)"' 
IREs from z=l to z=0.4, tire being the characteristic timescale. An efficient and episodic star formation is further supported by 
the luminosity-metallicity relation of z~ 0.7 emission line galaxies, which we find to be on average metal deficient by a factor 
of ~ 2 when compared to those of local spirals. We then examine how galaxy IREs can be related to the emergence at high 
redshift of the abundant population of galaxies with small size (but not with small mass), blue core and many irregularities. We 
show that recent merging and gas infall naturally explain both morphological changes and episodic star formation history in a 
hierarchical galaxy formation frame. 

We propose a simple scenario in which 75±25% of intermediate mass spirals have recently experienced their last major merger 
event, leading to a drastic reshaping of their bulges and disks during the last 8 Gyrs. It links in a simple manner distant and 
local galaxies, and gives account of the simultaneous decreases during that period, of the cosmic star formation density, of the 
merger rate, and of the number densities of LIRGs, compact and irregular galaxies, while the densities of ellipticals and large 
spirals are essentially unaffected. It predicts that 42, 22 and 36% of the IR (episodic) star formation density is related to major 
mergers, minor mergers and gas infall, respectively. 
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1. Introduction 

Understanding how and when galaxies and the stars in them 
were formed is still a considerable challenge for astrophysi- 
cists. A global picture of galaxy evolution during the last 8 
Gyrs has emerged from the results of all-inclusive galaxy sur- 
veys. The apparently unchanged number density of ellipticals 
and of large spirals at z~ 0.7 (Schade et al. 1999; Lilly et al. 
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1998; Conselice et al, 2004a) has lead to a picture for which 
most of the recent star formation occurred in small mass galax- 
ies (Cowie et al. 1996; Brinchman and Ellis, 2000, hereafter 
BE2000). This picture has been called "galactic downsizing" to 
illustrate its apparent contradiction to expectations from hier- 
archical models of galaxy formation. It is only consistent with 
a modest fraction of stars recently formed, since present-day 
stars are mostly locked in spheroids (bulge of spirals and ellip- 
ticals, Fukugita, Hogan & Peebles, 1998). 
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This is challenged by the evolution of the cosmic stellar 
mass density which is based on a detailed modelling of galaxy 
spectral energy distributions (Cole et al. 2001; Dickinson et 
al. 2003,; Fontana et al. 2003; Drory et al. 2004; Rudnick et 
al. 2003). Even if they suffer from various uncertainties (cos- 
mic variance, assumptions on stellar populations, metal abun- 
dances, extinction and initial mass function), all these studies 
predict that 30 to 50% of the mass locked in stars in present 
day galaxies actually condensed into stars at z < 1 . These es- 
timates agree with integrations of the universal star formation 
history, if it properly accounts for the Ught re-emitted by dust 
at IR wavelengths (Flores et al. 1999; Chary and Elbaz, 2001). 
Heavens et al. (2004, hereafter, H2004) have derived a peak of 
star formation at modest redshifts (z=0.5-0.6), mostly due to 
intermediate mass galaxies (3-30 IO'^'Mq) which have assem- 
bled/formed 60% of their stellar masses during the last 8 Gyrs 
(from integrations of their Figure 2). There is now a reasonable 
agreement that a significant part of the present-day stars (30% 
to 50%) was formed at recent epochs. 

Where are these relatively recent stars bom and how they 
have been formed ? We examine the major evolutionary trends 
since z=l from all inclusive surveys. Five to eight Gyrs ago, 
the galaxy population included 70% starbursts ( defined as 
Wo([OII])> 15A ), a fraction which has sharply decreased to 
17% in the local Universe (e.g. Hammer et al. 1997). Both UV 
(Lilly et al. 1996) and IR (Flores et al. 1999) comoving lumi- 
nosity densities have been declining rapidly since z = 1. 

1.1. Most of the star formation since z=1 is related to 

the IR luminosity density 

At IR wavelengths the comoving luminosity density was 70 
times larger at z = 1 than it is today (Elbaz et al. 1999; Elbaz 
and Cesarsky, 2003). At z = 1 it is mostly due to luminous IR 
galaxies (LIRGs, defined as Ljr > 10^ 'Lq following Sanders 
and Mirabel, 1996), which are mainly powered by star forma- 
tion (Fadda et al 2002). LIRGs form stars at very high rates 
(from 20 to several hundred of Mg/yr) and are responsible for 
the bulk of the cosmic infrared background (Elbaz et al. 2002). 
An early conclusion from studies of mid infrared emission and 
SEDs of galaxies is that from 50 to 66% of the CSFD at z ~ 0.7 
is related to far infrared emission (Flores et al. 1999), mostly 
due to bright mid-lR sources (LIRGs). The large impact of star 
formation revealed by the IR has been recently confirmed by 
Spitzer counts of 24 fan sources (Papovich et al. 2004). 

Distant LIRG morphologies encompass those of disks, ir- 
regulars and LCGs, with a noticeable fraction of them show- 
ing signs of merging or interactions (Zheng et al. 2004a). If 
sustained at the observed high rates, their star formation can 
double their mass in less than 1 Gyr (Franceschini et al. 2003; 
Zheng et al. 2004a). Liang et al. (2004b) have derived their ex- 
tinction, star formation rate and metal abundances from a de- 
tailed analysis of their high quality spectra from VLT. Distant 
LIRGs present a metal content less than half of that of local 
bright disks. Liang et al. (2004b) and Zheng et al. (2004a) have 
suggested that many distant disks are still forming a significant 
fraction of their stellar mass content. 



1.2. Can we relate the UV luminosity density decline 
to the evolution of the merging rate ? 

Peculiar galaxies are found to be the most rapidly evolving 
galaxy population and to be responsible for most of the de- 
cline of the UV luminosity density (Brinchman et al. 1998; 
Lilly et al. 1998, BE2000, ConseUce et al, 2004b). Besides 
showing unusual or irregular morphologies, peculiar galaxies 
generally display small half light radii (Lilly et al. 1998), while 
their stellar masses are comparable to those of large, more reg- 
ular galaxies (BE2000). It has been often suggested that the 
morphologies of peculiar galaxies are related to interactions 
and that merging could play a major role in their number den- 
sity decline. This is broadly consistent with the large expected 
amount of merging events during the last 8 Gyrs. From number 
counts of pairs up to z = 1, Le Fevre et al. (2000) estimated 
that each luminous (L*) galaxy experiences 1.1 to 1.8 x (7 10** 
yrs/Tmerger) major merging event (with a companion with half 
and one fourth of its luminosity, respectively). In this study the 
characteristic time over which the merger is completed, T^erger 
was assumed to range from 0.4 to 1 Gyr following expectations 
from simulations (Hemquist & Mihos, 1995). A recent study 
made in the near IR (Bundy et al. 2004) leads us to suspect a 
slower merger history, and accounting for the various selection 
eff'ects (colour & relative spatial resolution) would decrease the 
Le Fevre et al. estimates to half their original values. Merging 
appears still to be an important ingredient in galaxy evolution 
at z< 1, and this matches rather well with a recent stellar mass 
production. 

Several clues also favor a merger origin for many pecu- 
liar galaxies. In the CFRS, Hannmer et al. (2001) have stud- 
ied in detail the most compact ones {rso < 3.5 kpc), hereafter 
called luminous compact galaxies (LCGs) which show blue 
U-V colours but not necessarily blue B-K colours. LCGs are 
found to be dusty enshrouded starbursts (SFR averaged to 40 
Mg/jr) superimposed on a mix of metal-rich old stellar pop- 
ulation and intermediate age stars. Hammer et al (2001) sug- 
gested that LCGs could be remnants of mergers or interactions, 
and could be an important stage in the bulge formation in many 
spirals. Detailed studies of few local counterparts by Barton & 
Zee (2001) and by Osdin et al. (2001) have revealed tidally 
triggered central starbursts or complex velocity fields typical 
of completed mergers, which also explain their narrow emis- 
sion lines. There are also theoretical grounds which support 
that merger remnants could be observed in large numbers in the 
distant Universe. Although after a major merger, the dynami- 
cal relaxation time is short (few 10^ yrs, Barnes & Hernquist, 
1992), hydro-dynamical models (Tissera et al. 2002; see also 
Baugh et al. 1996) show that enhanced star formation can be 
sustained from several 10** yrs to ~ 2 10^ yrs after the collision, 
a period during which formation of stars in the surrounding 
disk can be accommodated. 

1.3. The present study 

Thus there are good reasons to suspect that merging probably 
has an important impact on the observed recent formation of 
stars. Because a significant part of the recent star formation is 
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related to dust-enshrouded galaxies, it calls for a panchromatic 
study of galaxy properties. This study is based on long-term 
follow-up observations of the CFRS by our team, using HST, 
ISO, VLA and VLT. It has the goal to hnk the morphological 
changes seen in optical wavelengths to the star formation re- 
vealed at IR wavelengths, and, more ambitiously, to be able 
to see an overall picture of galaxy evolution over the last 8 
Gyr. It targets 195 z > 0.4 galaxies, mostly from the sample 
of the CFRS which is probably the most widely studied sur- 
vey of galaxies for its statistical properties. CFRS is essentially 
complete up to z~ 1, encompassing all luminous (Mb < -20) 
galaxies (Hammer et al. 1997) with stellar masses from 3 to 30 
IO'^'Mq, those being responsible for the bulk of the star forma- 
tion history reported by H2004. There is a general consensus 
that intermediate mass galaxies dominate the stellar mass bud- 
get at z < 1 (80% of the total mass according to BE2000, 65% 
from integrations of H2004). 

In Section 2 we present the data and a comparison of the 
various estimates of star formation rate is given in Section 3. 
A soUd estimate of the stellar mass fraction formed in the past 
8 Gyrs is provided in Section 4, as well as new insights on the 
star formation histories of individual galaxies. Section 5 de- 
scribes the major morphological changes during this elapsed 
time, which are linked to the galaxy star formation history 
in Section 6, leading to a simple scenario which accounts for 
all changes in the intermediate mass galaxy population, in the 
framework of the hierarchical galaxy formation theory. ACDM 
cosmology (7/o=70 kms ' Mpc"^, Q.m = 0.3 and = 0.7) is 
adopted throughout the paper. 



2. The data obtained from our panchromatic study 
of distant galaxies 

This randomly selected sample of galaxies (for a detailed de- 
scription, see Liang et al, 2004b; Zheng et al, 2004a) was con- 
tructed in order to study each subclass of galaxies (LIRGs, 
LCGs, disk galaxies and starbursts). To study LIRGs two ad- 
ditional ISOCAM fields were included to avoid any contami- 
nation due to cosmological variance: selection and imagery on 
extra deep ISOCAM fields were done in the same way than for 
the CFRS survey (Iab < 22.5 galaxies). 
Table 1 summarizes the properties of 195 Iab ^ 22.5 galaxies 
with z > 0.4, on which 163, 14 and 23 are collected from the 
CFRS, the UDSR and the UDSF fields, respectively. Redshift 
measurements (and derived absolute magnitudes in the AB sys- 
tem) have been derived from CFRS or from VLT spectroscopy. 
The VLT follow-up was based on optical and near IR spec- 
troscopy' at resolution R> 1000, aiming at studying properties 
of z > 0.4 galaxies in the fields of the southern hemisphere 
(all fields but the CFRS14h). The new data gathered here al- 
low us to compare the physical properties of the most evolving 
galaxies (LIRGs and LCGs) to those of the rest of the galaxy 
population detected by the CFRS. In the following, conclusions 
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relative to the general galaxy population will be systematically 
derived from statistics based on the whole CFRS sample. 

The HST follow-up was based on numerous WFPC2 
archival data on the CFRS 03h and CFRS 14h fields on which 
we select galaxies with deep images in two colours (mainly 
F606W and F8 14W, see Zheng et al. 2004a for details). Central 
colours and compactness are derived from multiple color HST 
images (see Zheng et al. 2004a). 

Stellar masses of galaxies can be estimated from their near 
IR luminosities (Chariot, 1999; Bell and de Jong, 2000). Here 
we have chosen a conservative approach, which assumes that 
M/Lk depends on the rest-frame B-V color following the re- 
lationship derived by Bell et al (2003). We verified that an al- 
ternative use of the Bell & de Jong (2000) relationship would 
not much afffect our results and typical uncertainties are 0.1-0.2 
dex (Bell et al, 2003). Rest frame colors and absolute K band 
magnitudes were derived from observed V-I and I-K (and K 
magnitude), respectively. This leads to stellar mass values in 
the same range as those derived by BE2000 for the CFRS sam- 
ple: most stellar masses of z> 0.4 CFRS galaxies range from 
10.5 to 1 1.5 on a logarithmic scale. 

Four fields (CFRS 03h, 14h, UDSR and UDSF) were ob- 
served by deep ISOCAM pointings (Flores et al. 1999; 2004; 
Elbaz et al 2004) providing measurements of 15/im luminosi- 
ties which have been converted in IR luminosities (see Elbaz 
et al. 2002). Star formation rates based on 3 different luminosi- 
ties (IR, Ha, [Oil] 3727) assume the same IMF (see Hammer et 
al. 2003 and section 3). Ha luminosities are derived from VLT 
measurements (see Flores et al. 2004; Liang et al. 2004b) while 
some [Oil] 3727 luminosities are derived from CFRS spectra 
after a proper account for aperture corrections (see Liang et 
al. 2004b). 0/H gas abundances are calculated (see Liang et 
al. 2004b) from the ratio of [OII]3727 and [OIII]5007 fines 
to Hp after a proper account for extinction effects. A few 
values of 0/H are however directly derived from the single 
[OIII]5007/i?yS ratio (see Liang et al. 2004b). 

3. Towards more robust estimates of star 
formation in distant galaxies 

A prerequisite for estimating star formation from emission 
lines is to properly estimate the extinction. Ha luminosity is 
one of the best indicators of the instantaneous SFR. However, 
low resolution spectroscopy (R < 500) often produces mislead- 
ing results (Liang et al. 2004a). Only good S /N spectroscopy 
with moderate spectral resolution (R> 600) allows a proper es- 
timate of the extinction from the Ha/H/3 ratio after accounting 
for underlying stellar absorption. SFRs can be otherwise un- 
derestimated or overestimated by factors reaching 10, even if 
one accounts for an ad hoc extinction correction. These effects 
are prominent for a large fraction of evolved massive galaxies 
especially those experiencing successive bursts (A and F stars 
dominating their absorption spectra). 

ISOCAM observations have provided a unique opportunity 
to test the vaUdity of our SFR estimates. It has been shown 
(Elbaz et al. 2002) that the mid-IR and radio luminosities cor- 
relate well up to z = 1. Bolometric IR measurements derived 
from the mid-IR are validated, unless if both the radio-FTR and 
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Table 1. Sample of distant galaxies 
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N/D 








0.748 


0.43 


compact 


14.0422 


0.4210 


-20.08 


-22.29 


10.73 


N/D 




0.50 




2.078 


0.31 


e/sO 


14.0446 


0.6030 


-21.52 






89.03 








1.582 


0.44 


spiral 


14.0485 


0.6545 


-19.93 


-21.10 


10.18 


N/D 




5.20 




1.357 


0.22 


spiral 


14.0547 


1.1963 


-22.44 


-23.83 


11.28 


214.8 












merger 


14.0557 


0.8800 


-21.44 






N/D 








0.877 


0.59 


compact 



the MIR-FIR correlations have strongly changed at high red- 
shifts. Mid-IR photometry provides a unique tool to estimate 
properly the SFR of LIRGs (defined as Lir > 10^' Lq) up to 
z = 1.2, and of starbursts (Lir < lO'^ L©) up to z = 0.4. For 
a sample of 16 ISO galaxies at < z < 1, SFRs have been 
derived from the extinction-corrected Ha luminosities (Flores 
et al. 2004). These values agree within a factor 2 with mid-IR 
estimates (median value of S FRjr/S FR[j„-l.05). Using mod- 
erate spectral resolution (R= 1200) with a good S/N of 90 ISO 
galaxies (0.2 < z< 1), Liang et al. (2004b) have shown that the 



extinction can be properly derived from the HpjHj ratio, us- 
ing high S/N spectroscopy at moderate resolution at VLT. Such 
measurements can be performed for all galaxies up to z = 1, 
and Liang et al (2004b) find that the derived SFR from extinc- 
tion corrected Hj5 line is consistent with S FRir, with their ratio 
ranging from 0.6 to 1.7. 

Can we also use [Oil] 3727 luminosity to derive SFRs of 
high redshift galaxies using spectrographs in the visible range 
? Figure 1 shows that it systematically underestimates the SFR 
by large factors, reaching median values of 5 to 22 for star- 
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CFRS 


z 


Mfi 




log(M/Mo) 


IR 


SFR^ 

^ //a 


fOU] 


12+log(0/H) 


c.c." 


Comp*^ 


Morph'' 


14.0580 


0,7440 


-21 07 






N/D 








1.849 


0.28 




1 4 OfiOO 


1.0385 


-21.75 






636.2 




27.75 








lITCg ulsX 


14.0651 


0.5480 


-19.44 


-21.75 


10.51 


N/D 








1.674 


0.53 


e/sO 


14.0663 


0.7434 


-21.40 






187.5 




2.89 




1.580 


0.37 


spiral 




0.8110 


-19.89 


-20.85 


10.07 


159.5 














14.0697 


0.8271 


-20.71 






N/D 








0.914 


0.69 


compa.ct 


14.0700 


0.6526 


-21.29 






N/D 








1.905 


0.65 


e/sO 


14.0711 


1.1 180 


-22.09 


-23.30 


1 1 .07 


N/D 








1.517 


0.41 


irrcguliir 


14.0725 


0.5820 


-19.54 


-21.70 


10.47 


41.53 




3.12 




1.543 


0.26 


irrcgLila.r 


14.0727 


0.4638 


-20.33 






N/D 




7.88 




1.804 


0.35 




14,0746 


0.6750 


-20.40 


-22.64 


10.88 


N/D 








1.968 


0.50 


e/so' 


14.0814 


0.9995 


-20.50 


-22.87 


10.90 


N/D 








1.636 


0.17 


iirpCTiilat" 


14.0820 


0.9800 


-21.53 


-24.08 


11.45 


N/D 








2.300 


0.43 


e/sO 


14.0846 


0.9900 


-21.37 


-22.99 


10.96 


N/D 








1.353 


0.25 




14.0854 


0.9920 


-21.66 


-23.83 


11.34 


N/D 








2.021 


0.45 


e/sO 


14.0899 


0.8750 


-21.07 






N/D 




4.69 




0.438 


0.34 




14.0913 


1 .0058 


-21.58 






N/D 








1.946 


0.58 


e/so' 


14.0937 


1.0097 


-21.60 






N/D 








2.248 


0.47 


e/sO 


14.0962 


0.7617 


-21.30 


-23.02 


11.02 


N/D 




4.30 




1.900 


0.52 


spiral 


14.0964 


0.4351 


-19.32 






N/D 








1.232 


0.27 




14.0972 


0.6809 


-21.10 


-22.02 


10.55 


N/D 




27.06 




0.880 


0.54 




14.0998 


0.4300 


-20.11 


-22.23 


10.70 


44.90 








1.795 


0.42 


itTPCTiilni" 


14.1008 


0.4362 


-19.82 






N/D 








2.475 


0.43 




14. 1012 


0.4815 


-19.85 


-20.92 


10.13 


N/D 




2.36 




1.248 


0.42 


conipact 


14.1028 
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-21.75 


-23.67 


11.28 


N/D 




14.48 




2.060 


0.39 


e/sO 


14.1037 


0.5489 


-20.31 


-21.50 


10.37 


N/D 




2.45 




1.381 


0.17 




14.1041 


0.4351 


-19.50 






N/D 












irregular 


14.1042 


0.8916 


-21.42 


-23.30 


11.12 


161.6 




2.33 




1.773 


0.42 




14.1043 


0.6516 


-21.60 


-23.64 


11.26 


N/D 








2.035 


0.42 




14.1079 


0.9130 


-21.01 


-20.95 


10.15 


N/D 




8.76 




1.633 


0.33 


irregular 


14.1087 


0.6595 


-20.21 


-21.08 


10.17 


N/D 




5.34 




1.061 


0.44 


e/sO 


14.1129 


0.8443 


-21.19 


-22.61 


10.81 


21 1.9 












merger 


14.1136 


0.6404 


-20.16 


-19.12 


9.38 


N/D 




10.73 




0.723 


0.43 


compact 


14.1139 


0.6600 


-21.60 


-23.45 


11.15 


203.2 




7.70 








merger 


14.1146 


0.7437 


-20.77 






N/D 




8.83 




1 .450 


0.51 


coniptict 


14.1157 


1.0106 


-22.49 






N/D 












merger 


14.1164 


0.6773 


-20.54 






N/D 








1.066 


0.30 


itTPCTiilat" 


14.1179 


0.4345 


-19.13 






N/D 








1.953 


0.50 


e/sO 


14.1189 


0.7526 


-20.42 


-21.57 


10.39 


N/D 




5.33 




1.143 


0.35 




14.1190 


0.7544 


-21.38 






159.1 




3.06 




1.583 


0.35 


spird 


14.1213 


0.7631 


-21.03 






N/D 














14.1223 


0.6936 


-20.42 






N/D 








1.078 


0.32 




14.1232 


0.7660 


-20.22 


-22.73 


10.85 


N/D 












merger 


14.1246 


0.6486 


-19.76 






N/D 








2.140 


0.27 


irregular 


14.1277 


0.8100 


-21.33 






N/D 




6.91 




1.917 


0.52 




14.1305 


0.8069 


-20.59 






N/D 












m.erger 


14.1311 


SOfiS 


-22.03 






N/D 








2.022 


0.52 


e/sO 


14.1326 


0.8111 


-20.74 






N/D 








1.200 


0.20 


irregular 


14.1329 


0.3750 


-20.63 






22.43 




3.48 










14.1350 


1 .0054 


-21.29 






N/D 








1.631 


0.18 


irregular 


14.1355 


0.4801 


-19.74 


-20.63 


10.01 


N/D 




3.94 




1.412 


0.26 


.p 


14.1356 


0.8307 


-20.44 






N/D 




7.61 




1.739 


0.42 


compaci 


14.1386 


0.74 1 


-21.22 


-22.72 


10.84 


132.32 




8.50 










14.1395 


0.5301 


-19.85 






N/D 




4.47 




1.070 


0.27 


spiral 


14.1412 


0.9122 


-21.03 






N/D 








1.729 


0.25 


Spiral 


14.1415 


0.7486 


-21.15 






N/D 








2.014 


0.48 


compact 


14.1427 


0.5379 


-19.77 






N/D 








1.459 


0.42 


spiral 


14.1444 


0.7420 


-19.89 


-22.57 


10.82 


139.57 




2.91 










14.1464 


0.4620 


-19.62 






N/D 




0.69 




2.014 


0.56 


e/sO 


14.1477 


0.8191 


-21.02 






N/D 








2.002 


0.52 


e/sO 


14.1496 


0.9163 


-21.06 






N/D 




14.01 




1.097 


0.55 


compact 


14.1501 


1.0018 


-21.41 






N/D 












irregular 


14.1506 


1.0126 


-21.51 






N/D 








1.942 


0.59 


e/sO 


14.1524 


0.4297 


-21.16 






N/D 








1.363 


0.31 


spiral 


14.1525 


0.7480 


-20.62 


-22.00 


10.55 


250.90 




13.10 











bursts and for LIRGs, respectively. The fact that UV fluxes 

provide even lower SFR values than [Oil] 3727 fluxes might 
be related to the expected increase of the [OII]3727/7/yS ratio 
with decreasing metal abundance (high z systems being ex- 
pected to be less metal abundant than local ones). An important 
word of caution is however required here, since using the UV 
slope, Meurer, Heckman & Calzetti (1999) were able to derive 
extinction corrections after comparison to modelling. By com- 
paring these corrections to IR fluxes of z~ 0.9 ISO galaxies 
from Flores et al. (1999), Adelberger & Steidel (2000) claimed 



that their method appears to work within a factor 3 (see also 
Adelberger, 2001). On the other hand. Bell (2002) argues that 
the rest-frame UV and optical data provides estimates of the 
star formation rates which are subject to large systematic un- 
certainties because of dust, which cannot be reliably corrected 
for using only UV/optical diagnostics. 

In summary, the star formation for strong star forming sys- 
tems is severely underestimated by using UV tracers alone, 
which is not surprising given the expected effects of dust ex- 
tinction. It is hardly a new result since the discovery of LIRGs 
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CFRS 


z 


Mb 


Mk 
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0.5370 
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0.39 
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0.6243 


-19.93 


-21.72 
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0.5419 
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-21.37 
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0.4070 


-18.83 


-20.46 


9.99 
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0.9252 
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-22.37 
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31.66 


6.38 


8.58 








22.0828 


0.4070 


-18.75 


-20.53 


10.01 
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1.25 
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0.4714 


-19.92 
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0.83 


8.17 
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0.5369 


-19.76 


-21.64 


10.40 
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5.65 


8.46 
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0.8796 


-20.49 
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8.55 










0.4501 
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15 
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17 
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20 
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26.64 
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21 
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23 
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26 
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3.17 
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27 


0.6429 


-20.67 
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30 
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3.41 
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14.98 
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10.19 
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02 
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4.14 


8.56 








03 


0.5532 


-20.52 






1 1.14 






8.82 








04 


0.9620 


-21.55 






97.17 




3.17 










06 


0.6928 


-21.49 






15.97 




3.46 










07 


0.7014 


-20.70 






16.77 


11.43 


2.86 


8.81 








08 


0.7075 


-21.03* 






14.01 




2.19 










12 


0.7388 


-20.95 
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5.99 
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-20 07 
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7 70 
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16 
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17 


0.8100 
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38.18 
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18 


0.4620 


-20.60* 
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19 


0.5476 


-21.94* 






80.86 


60.94 


4.71 


8.92 








21 


0.6980 


-21.35* 






16.20 


14.40 




8.74 








24 


1.1590 


-23.31* 






138.9 




57.02 










25 


0.8094 


-20.42* 






12.52 




11.46 










26a 


0.7023 


-21.60* 






16.31 


10.45 


12.68 


8.36 








26b 


0.7027 


-20.62* 






18.38 


27.11 




8.62 








28 


0.6621 


-21.84* 






42.96 














29 


0.6619 


-20.45* 






51.34 


23.34 


2.44 


8.43 








31 


0.6868 


-20.57* 






21.71 




1.56 


8.82 








32 


0.7268 


-19.73* 






32.28 




1.25 











— Not available 
N/D Not detected 

" Star formation rate in Mo/yr. 
c,d & e Central color, compactness and morphology from HST images 
* Mb estimated from observed R band and from the observed spectroscopic slope 



& ULIRGs by IRAS. We will see in the following that, unlike 
the situation in the local Universe, such an underestimation has 
profound implications for our understanding of the evolution in 
the distant Universe. 



4. How and when did galaxies form the bulk of 
their stars and metals ? 

4.1. Luminous IR galaxies: a rare phenomenon locally, 
a common one atz> 0.4 

LIRGs are rare objects in the local Universe (Sanders and 
Mirabel, 1996). Soifer et al. (1986) estimated that 0.5% of 
galaxies with visible luminosities larger than IO^^Lq, have 
spectral energy distribution typical of LIRGs (90% of their 
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Fig. 1. Comparison of SFRs calculated using the formalism of Kennicutt (1998). (Left panel): [OII]/l3727 estimates are compared 
to IR estimates for a sample of 70 galaxies from the CFRS and Marano fields (see Liang et al. 2004b). Upper left panel provides 
the S FRiji/S FR[oii] ratio. Median values of the ratio are 5 and 22 for starbursts (SFR < 20 Mojr"') and LIRGs respectively. 
(Right panel): same comparison for 2800A UV luminosity estimated for 61 CFRS galaxies. Median values of the S FRjr/S FRuv 
ratio are 13 and 36 for starbursts and LIRGs, respectively. 



luminosity emitted in IR). How this compares to the distant 
Universe ? 

A total of 113 galaxies have been detected at 15yum by 
ISOCAM (Flores et al. 1999, 2004) in the CFRS 14h and 03h 
fields (2 X 100 arcmin^), and 70% of them are LIRGs at z > 0.4, 
all with Mb < -20. Compared with the 311 z > 0.4 CFRS 
galaxies occupying a solid angle of 112 arcmin^ (Hammer et 
al. 1997), it appears that 15% of the Mb < -20 galaxies at 
z > 0.4 are indeed LIRGs. A similar calculation can be made by 
considering the subsample of 184 galaxies possessing K band 
photometry, providing a similar fraction (15%) of LIRGs in M 
> 5 10"^ Mq galaxies at z > 0.4. This fraction is remarkably 
stable to varying the lower mass limit from 3 10'" Mq to 10 
10'" Mq. The marginal incompleteness of the CFRS (85% of 
redshift identified) can slightly decrease the LIRG fraction and 
it is mostly related to the difficulty in deriving a redshift for 
distant red galaxies near the CFRS Iab limit (Crampton et al. 
1995; Hammer et al. 1997). On the other hand, the fraction of 
LIRGs is underestimated because the depth of the ISO observa- 
tions of the CFRS fields (complete at S(15yum) > 300 yuJy) does 
not ensure that all LIRGs near the cut-off limits (L/« = 1O"L0) 
have been detected, especially at the high-z end (see the dis- 
cussion by Liang et al. 2004b). A gross consistency check is 
provided by comparing the mid-IR counts from Elbaz et al. 
(1999) to those from optical, knowing that the redshift distri- 
bution of mid-IR faint sources is roughly comparable to that of 
the CFRS. The number density of lOOyuJy sources is 31% that 
of Iab < 22.5, and the Franceschini et al. (2003)' study of such 
sources in the HDF South find that 11/21 of them are indeed 



LIRGs with Iab < 22.5. This is independent support that at z> 

0. 4, 15% of the CFRS galaxies are LIRGs. 

4.2. Stellar mass formed during the time between 
z ^ 1 andz ^ 0.4 

The ratio of stellar mass to SFR provides the time-scale to dou- 
ble the stellar mass content, assuming a constant SFR. Figure 

2 relates the doubling timescale to the stellar masses of LIRGs 
and LCGs. Re-evaluating the SFRs of galaxies from IR data 
leads to doubling time-scales from 0.1 to 1.1 Gyr for LIRGs, 

1. e. much shorter than those derived in earlier work (BE2000). 
Here we investigate how the new stars formed in LIRGs can 
impact the stellar mass of the galaxy population in the interme- 
diate redshift range. 

The intermediate stellar mass galaxies (10.5 < 
log(MIMQ) < 11.5) include 15% of LIRGs with a me- 
dian doubling time-scale of 0.8 Gyr (Figure 2), in excellent 
agreement with the median value from Franceschini et al. 
(2003, see their Figure 10) for S(15//m) > 100 yuJy galaxies 
within the same mass range. The newly formed stellar mass 
in LIRGs during the elapsed time of 3.3 Gyr (from z = 1 to 
z = 0.4) correspond to 0.15 x 3.3/0.8 = 62% of the z = 1 total 
mass of intermediate mass galaxies. Flores et al. (1999) found 
that by including the IR contribution, the CSFD is between 2 to 

3 times larger than that estimated from UV. Including all kind 
of star formation (IR & UV emittors), it is expected that the 
recent star formation produces from 70% to 100% of the total 
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Fig. 2. Typical timescale for forming stellar mass in z > 0.4 LIRGs (full dots) and LCGs (filled squares, no LIRGs included) 
related to their observed stellar masses. SFRs have been estimated from extinction-coiTected Balmer emission lines (LCGs) or 
IR luminosities (LIRGs). The dotted line box indicates the position occupied by the galaxies when their SFR is derived from 
[OII]3727 line observations, following the adopted calibration from Brinchmann and Ellis (2000). Similar estimates have been 
made for this sample for LIRGs (empty dots) and for compact galaxies (empty squares). This shows that an important part of the 
star formation can be hidden using UV estimates of the SFR. The horizontal dashed line shows the median doubling time-scale 
for the population of intermediate mass LIRGs. There is a hint that LIRG masses could have been overestimated by a factor of ~ 
2 (see section 7); this would move the full dots towards the upper-left, leading to a median characteristic time of 0.4 Gyr, while 
the total mass formed would be unchanged (SFR still stands). Accurate determinations of SFR/M for LIRGs will be provided by 
Marcillac et al. (2004, in preparation) 



mass of intermediate mass galaxies at z = 1 (or alternatively 
40 to 50% of the mass at z = 0). 

The star formation in LIRGs provides a substantial stel- 
lar mass production in intermediate mass galaxies, because of 
their large number density at high redshifts, in constrast with 
the situation in the local Universe. The extremely tiny fraction 
of recently formed mass by BE2000 is due to their use of the 
[Oil] 3727 line as a SFR tracer, and they indeed commented 
that this choice could lead to serious underestimates. Our re- 
sult is however in good agreement with methods independent 



of SFR estimates from UV tracers, such as those based on fos- 
sil records of present-day intermediate mass galaxies (H2004). 
Integrations from z - 1 to z = 0.4 of their SFR history lead to 
a stellar mass increase by ~110% of their total mass at z = L 
A non-negligible mass accretion from satellites (see Bundy et 
al. 2004) could also contribute to the stellar mass assembly in 
intermediate mass galaxies. 
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Fig. 3. Median combined VLT/FORS (R600+I600) spectra of 38 distant LIRGs, 21 LCGs (LIRGs have been excluded here) and 
26 normal starbursts showing their absorption lines. A 3cr clipping method has been used to reject residual unrelated features 
(mostly from sky substraction). All spectra invariably show a mixture of intermediate age stellar population revealed by their 
strong Balmer absorption system with an older stellar population revealed by metallic absorption lines. 



4.3. Successive infra-red episodes (IREs): a plausible 
star formation history for individual galaxies 

How can we interpret such a large number density of LIRGs 
in the distant Universe ? Let us examine an extreme hypothesis 
that LIRGs are a specific galaxy population, and are actively 
forming stars at very high rate during the 3.3 Gyr from z = 1 
to z = 0.4. During that time these galaxies would multiply their 
masses by 2 x (3.3/0.8) = 8.2, producing galaxies in the mass 
range 11.4 < logiM/Mo) < 12.4, which is dominated by mas- 
sive ellipticals (Nakamura et al. 2004). Besides the theoretical 
challenge in interpreting a sustained star formation at such high 
rates (~ SOM^/yr), recent formation of most of their stars in a 
noticeable fraction of massive ellipticals in the field has never 
been detected by studies of their local properties. It is generally 
argued that they are formed at a much higher redshifts, possibly 
by monolithic collapses (Renzini, 1999). 



The only plausible alternative is that LIRGs are related to 
the population of intermediate mass galaxies, probably through 
an episodic star formation history. Hierarchical models in the 
ACDM cosmogony predict that galaxies formed through a 
number of mergers leading to star formation histories with mul- 
tiple short-living peaks (Abadi et al. 2003). A significant frac- 
tion of LIRGs are indeed associated with mergers at different 
stages after careful examinations of their morphologies (Flores 
et al. 1999; Aussel et al. 1999; Zheng et al. 2004a). Moreover, 
distant LIRGs display optical properties astonishingly similar 
to those of other distant field galaxies. Optical properties are 
unable to properly distinguish LIRGs from the more quies- 
cent field population. As for other distant galaxies, their spec- 
tra (Figure 3) invariably show the combination of an old stellar 
population, with strong Balmer absorption systems (A-hF stars) 
and emission lines. 
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Fig. 4. Mb versus OIH for CFRS galaxies (16 LIRGs, 12 LCGs and 19 starbursts at z> 0.4). Most low-z sources (z < 0.4) lie 
near the local spiral sequence, in contrast with z > 0.4 galaxies which on average, are metal deficient by a factor of 2 (see text). 
Data for local spirals and for z< 0.4 galaxies are described in Liang et al (2004b). 



Therefore, the most plausible explanation of the high frac- 
tion of LIRGs in the distant Universe is that they correspond 
to episodes of violent star formation that are superimposed 
on relatively low levels of (secular?) star formation activity 
in more quiescent galaxies. It is generally believed that typi- 
cal timescales for strongly enhanced star formation in LIRGs 
are ~ O.lGyr (Franceschini et al. 2003), and this is consis- 
tent with peaks of star formation occuring during merging 
events described by hydrodynamical models (Tissera et al. 
2002; Scannapieco & Tissera, 2003). 

If we assume that all intermediate mass galaxies could po- 
tentially experience an infrared episode (IRE), during a char- 
acteristic time Tire, then the number of episodes per galaxy 
is niRE - 0.15x (3.3G}'r/TiRE) - 5(TiRE/0.1Gyr)"'. A lower 
limit on tire is provided by the lifetime of massive stars and 
dust formation (typically 0.01 Gyr). Enhanced star formation 
during merging episodes could have durations of several tenths 
of Gyr (see a description in Tissera et al. 2002). An impor- 
tant support for tjre ~0.1 Gyr is provided by Marcillac et al. 
(2004, in preparation) who can reproduce the Balmer absorp- 
tion line H8 and the 4000A break of LIRGs using a large set of 
star formation histories from Monte Carlo simulations based on 



Bruzual and Chariot (2004) models. Several infrared episodes 
per galaxy from z = 1 and z - 0.4 can easily account for the 
large number of z > 0.4 CFRS galaxies which show an impor- 
tant population of intermediate age stars (Hammer et al. 1997). 

4.4. The luminosity-metal abundance relation 

The gas phase metal content of 16 distant LIRGs and 12 LCGs 
was derived from VLT spectra (Hammer et al. 2001; Liang et 
al. 2004b; Gruel, 2002). The same methodology was applied 
(see Liang et al. 2004b for a detailed description) to evaluate 
the OjH ratio from emission lines of 19 additional z > 0.4 
normal starbursts, which form fewer stars than LIRGs and are 
more extended than LCGs. Figure 4 shows that all exhibit a 
very broad distribution of their gas abundances. There is a hint 
that LIRGs show slightly higher metal abundances than the rest 
of the sample, which is not unexpected given the high produc- 
tion of metals in such strong star formation sites. Such an ex- 
pectation might however be blurred by the 0.08 dex typical er- 
ror bars (Liang et al. 2004b). 

On average z > 0.4 galaxies show metal abundances two 
times lower than those of the local spiral sequence (Figure 
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Table 2. Morphological classification statistics for 36 LIRGs 
and 75 non LIRGs (z > 0.4), based on state of the art automatic 
classification schemes (Zheng et al. 2004a; also Zheng et al. 
2004b); for comparison, the last column shows the fractions 
derived from Nakamura et al (2004) for the same galaxy mass 
range. 



Type 


z>0.4 


z > 0.4 


z>0.4 


local 




LIRGs 


non-LIRGs 


galaxies 


galaxies 


E/SO 


0% 


27% 


23% 


27% 


Spiral 


36% 


45% 


43% 


70% 


LCG 


25% 


17% 


19% 


< 2% 


Irregular 


22% 


7% 


9% 


3% 


Major merger 


17% 


4% 


6% 


<2% 



4). If, as we propose later (see section 6), most intermediate 
mass galaxies at z > 0.4 are connected to present-day spi- 
rals, this provides independent support that they (at least those 
with emission lines) have been actively forming about half of 
their present-day stellar and metal content from z ~ 1 to 0.4. 
All trends in Figure 4 can be easily accounted for by episodic 
star formation (displayed by LIRGs): the broad range of metal 
abundance of LIRGs takes its origin from that of the under- 
lying population of emission line galaxies. The dashed line in 
Figure 4 represents a simple infall model, assuming a 1 Gyr 
infall time and a total mass of 10'^ Mq (see Kobulnicky et al. 
2003). A more complex evolutionary track probably applies to 
distant emission line galaxies, as illustrated in the top left of the 
figure. After a LIRG episode the luminosity decreases because 
of the disappearance of short-lived massive blue stars. During a 
LIRG episode the ISM is enriched and the luminosity increases 
due to the burst. The ampUtude of both fading and burst bright- 
ening is however attenuated by the dust produced by massive 
stars. Dust extinction in LIRGs is severe enough that the UV 
luminosity of LIRGs is, on average, only ~2 times greater than 
that of starbursts (see Figure 1), which is consistent with their 
relative M^ distributions. 

Keck and VLT studies (Kobulnicky et al. 2003; Maier, 
Meisenheimer and Hippelein, 2004) find a similar metal de- 
ficiency for z > 0.4 galaxies as faint as M^ = -18. Taken 
together, these results indicate that, as an ensemble, emission 
line galaxies were poorer in metals at z ~ 0.7 than present day 
spirals by a factor of order two. 

5. Morphological changes 

5.1. Galaxy classification 

Distant galaxy morphologies can no longer be described by 
the Hubble classification sequence: many of them show irregu- 
lar, complex or compact morphologies (Brinchman et al. 1998; 
Lilly et al. 1998; Guzman et al. 1997). Table 2 shows the result 
of a morphological classification performed for 1 1 1 z > 0.4 
galaxies (see also Table 1) for which we have been able to 
derive their colour maps (see Zheng et al. 2004a for details 
on the procedure). Here we did not try to disentangle ellipti- 
cals from lenticulars; we have counted spirals as a single class 
and assessed another class (compact) for all objects barely re- 



solved by the HST. From their visual examinations, Brinchman 
et al (1998) generally classified our compact galaxies as pecu- 
liar (irregular & tadpole), while van den Bergh (2001) clas- 
sified most of them as spirals. Such a disagreement between 
the two studies of the same objects leads us to suspect that the 
classification of compact galaxies is not secure at the WFPC2 
spatial resolution. Our very drastic classification scheme based 
on deep images with two colours has the potential advantage 
of considerably limiting the uncertainties related to cosmologi- 
cal dimming, spatial resolution and morphological k-correction 
(see Abraham and van den Bergh, 2001). Apart from compact 
galaxies, our final classification is in excellent agreement with 
former ones (Brinchman et al. 1998; van den Bergh et al. 2001). 

Following BE2000, let us assume that the CFRS sample 
is somewhat representative of intermediate mass galaxies be- 
tween z = 0.4 and z = 1, with the morphological type distri- 
bution given in Table 2. Morphological determined luminos- 
ity functions (MDLF) in the local Universe are available from 
Nakamura et al. (2004), and one can approximate their mass 
function using the stellar mass/light ratio in the Sloan r-band 
from Glazebrook et al (2003). It provides for local intermedi- 
ate mass galaxies, a number fraction of E/SO and of early type 
spirals of 27% and 53%, respectively, the rest being mostly late 
type spirals. It is encouraging to see that, for early type galax- 
ies, the Nakamura et al (2004) MDLF does not differ substan- 
tially from that from Marzke et al (1998). It is important to no- 
tice that these fractions are barely affected (by a few percent) 
if we extend the lower mass Umit towards a higher value. The 
most remarkable trend is the small fraction, in the Nakamura 
et al sample, of peculiar galaxies such as irregulars, LCGs or 
mergers. An important goal of this paper is to approximately re- 
late the distribution of morphological types of Table 2 to those 
found in the local Universe. The most significant changes are, 
by decreasing order of significance, the complete vanishing of 
LCGs in the local Universe (factor =10 according to Jangren et 
al. 2004; see also Garland et al. 2003) and the number decrease 
of mergers and Irr (factor ~ 4). The fraction of E/SO appears 
reasonably constant with the epoch in agreement with Smith 
et al (2004) and the very modest increase of the E/SO fraction 
(also noticed by BE2000) is probably mainly due to the CFRS 
lack of completeness in measuring redshifts of red galaxies (see 
section 4.2). 

5.2. Strong star formation in central cores ofz> 0.4 
galaxies 

An important support for the idea of a recent evolution is also 
provided by the dramatic change in the galaxy central colours. 
Present-day galaxies in the Hubble sequence exhibit a red bulge 
superimposed on a blue disk. This is not longer the case for 
z > 0.4 galaxies. Figure 5 shows the distribution of their 
colours in the central 1 kpc region. A full methodological de- 
scription can be found in Zheng et al (2004a), and includes the 
use of V and I images superimposed with an accuracy of 0.015 
arcsec, a prerequisite for studying the colour of galaxies with 
complex morphologies. Almost all distant galaxies (including 
all LIRGs, Zheng et al. 2004a) have much bluer central colours 
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than local bulges, with the noticeable exception of distant el- 
liptical galaxies, which show similar colours to local ellipti- 
cals. We interpret this as evidence for an active star formation 
in the 1 kpc central region of most distant spirals, which can 
be either associated with star formation in their bulges, or with 
an absence of bulge at high redshift, or a combination of both 
effects. This supports a relatively recent formation of bulges 
(see Zheng et al., 2004b) in many present-day spirals which 
are mostly early type. 

6. Towards a scenario of a recent bulge and disk 
formation in intermediate mass spirals 

6.1. The need to up-date the current scenario of 
galaxy formation af z < 1 

Recent observational progress has lead to identification of more 
than a few rather robust features of the evolution of intermedi- 
ate mass galaxies. Up to now, no unanimity has been reached in 
interpreting them within a single picture. A "galactic downsiz- 
ing" picture (Cowie et al, 1996; BE2000) requires significant 
improvements, since apart from a recent and active formation 
of dwarves, it assumes that most galaxies have completed the 
bulk of their stellar mass at z = 1. Its major weaknesses are 
that: 

- the stellar mass formed since z = 1 (from 30 to 50% of the 
total) is unlikely produced in dwarf galaxies, because these 
objects contribute only a little to the local stellar mass or 
metals (Fukugita et al. 1995; H2004). 

- it neglects the role of dust-enshrouded star formation, the 
rate density of which has been found to be higher (Flores 
et al. 1999) than that derived from rest-frame UV light, and 
which is also required for matching the stellar mass evolu- 
tion after integrating the CSFD (see Dickinson et al. 2003). 

- it is not consistent with the reported deficiency of metal 
in the gas phase of intermediate mass galaxies at z > 0.4 
(Figure 4). 

- it predicts that intermediate mass spirals and ellipticals 
were mostly evolving passively since z = 1, while from 
pair counts, a significant fraction of L* galaxies should ex- 
perience major merger events between z = 1 and z = 
(section 1.2). 

- it does not easily account for the large number of pecu- 
liar galaxies (LCGs, irregulars & mergers) which repre- 
sent more than a third of the intermediate mass galaxies at 
z > 0.4, and display the most remarkable changes in galaxy 
morphology during the last 8 Gyr (Table 2). 

- it is qualitatively at odds with the predominance of blue 
cores in the distant galaxies which are not ellipticals. 

Of course many of the points above have been widely 
discussed in the literature, and some could have been partly 
solved in the framework of the "galactic downsizing" picture. 
Here we argue that the simultaneous changes in galaxy 
morphologies, central colours, stellar mass production and 
metal abundance properties are not fortuitous, and we propose 
to update the current scenario of galaxy formation during the 



last 8 Gyrs. To do this, we need to puzzle out the impact of 
the star formation history on the galaxy morphologies and to 
better understand how distant giilaxies are connected to the 
present-day Hubble sequence. 



6.2. Where and how have most stars formed since 
z=1 ? 

Almost all the recent star formation has occurred in spirals, 
LCGs, irregulars and major mergers, while the contribution of 
E/SO is negligible, both from UV (Lilly et al., 1995; Brinchman 
et al., 1999; Menanteau, Abraham & Ellis, 2001; Wolf et al., 
2004) and from IR lights (Zheng et al., 2004a; see also Table 2). 
In which present-day galaxy type has star formation occurred 
? Some of the major mergers (those having consumed all their 
gas) and some of the LCGs possibly gave rise to a recent forma- 
tion of a few intermediate mass ellipticals. However if the fate 
of all LCGs (or all LIRGs) were to become ellipticals, the den- 
sity of ellipticals today would be much higher than observed 
(see Table 2). Indeed many LIRGs show disk morphologies, 
many LCGs show faint surrounding disks, and so we assume 
that most of the recent star formation has occurred in progeni- 
tors of the numerous population of present-day spirals. 

The episodic star formation seen from IR data (see section 
4) is sufficient to explain most of the stellar mass formed since 
z = 1 , and is highly consistent with predictions from hierarchical 
galaxy formation models. In our view, LIRGs display the most 
active changes in galaxy morphology, while other galaxies are 
rather witnessing a more secular galaxy evolution. Zheng et al. 
(2004a) show that LIRGs present a well-defined sequence in 
the compactness versus central colour plane, from very blue 
compact cores (the LCGs) to grand design spirals with central 
colour consistent with present-day bulges. This observation is 
taken as the starting point to develop our galaxy formation sce- 
nario. Figure 6 illustrates an evolutionary sequence which as- 
sumes that several present-day spirals have been shaped dur- 
ing the recent formation of stars and metals in LIRG episodes. 
Indeed this sequence ends by an inside out formation of the 
disk, occurring from the outskirts of an already existing bulge. 
What we are proposing here is that, provided that new gas is 
available to the galaxy (in an internal reservoir or accreted from 
the vicinity), it re-emerges as a spiral. 

In their detailed simulations of the formation of a single spi- 
ral galaxy, Abadi et al. (2003) have assumed a relatively high 
redshift (z = 1.5) for the last major merger event, which pro- 
duces a short and intense starburst episode as well as a com- 
plete re-shaping of the galaxy bulge and disk. Our scenario 
is consistent with a more recent occurrence of the last major 
merger in most intermediate mass spirals. Galaxy mergers are 
very efficient in activating star formation as shown by hydro- 
dynamical models, and the (episodic) star formation history of 
our sequence can follow what it is displayed in Figure 1 of 
Tissera et al (2002). Major mergers also induce a shock-heated 
gas wind and the formation of a remnant gas halo, which could 
have a cooling time of more than one Gyr and could provide an 
enriched medium with which a gas disk could subsequently be 
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Fig. 5. Histogram of the 1 kpc inner central color, following Zheng et al. (2004a), for z > 0.4 galaxies including ellipticals 
(shaded area). Central colours have been calculated at z = 0.75, to minimize k-corrections errors (Hammer et al. 2001). The 
typical colour of a local bulge at z = 0.75 (assuming a passive evolution) is indicated. Morphologically selected ellipticals show 
no evidence for evolution while the core of most other galaxies are much bluer than expectations from a bulge passive evolution. 



formed (Cox et al, 2004). In the following we will test to what 
extent a scenario of "spiral rebuilding" is able to reproduce the 
evolutionary trends of intermediate mass galaxies. 



6.3. Ingredients for an updated scenario of "spiral 
rebuilding": consistency checks and predictions 

The sequence displayed in Figure 6 means that major merg- 
ers/compact/spirals & irregulars correspond to the various 
phases undergone by spiral progenitors from z-l to z=0.4. 
Let us also assume that the characteristic times correspond- 
ing to each phase can be approximated by the number frac- 
tion displayed in Table 2 (column 3). Table 3 shows these frac- 
tions rescaled after excluding E/SO, as well as the expected 
characteristic times if all (or alternatively 75% of) present- 
day spirals were formed at z<l following our evolutionary se- 



quence. The derived characteristic times match rather well with 
theoretical expectations based on simulations of major merg- 
ers. Our evolutionary sequence is then consistent with the ob- 
served fractions of the different galaxy types. Encouraged by 
this, we also consider the constraint from pair statistics (Le 
Fevre et al., 2000, revisited by Bundy et al., 2004, see section 
1.2). Assuming that 75+25% of galaxies experienced a major 
merger since z=l provides a particularily good match of the 
timescales in Table 3. 

Table 4 presents a score card showing different scenarios, 
including "galactic downsizing" and our scenario ("spiral re- 
building"), and compares their relative merits in reporting eight 
features of galaxy evolution since z=l. By construction our 
scenario is consistent with the CSFD evolution derived from 
UV and IR lights, with the IR background, and with the merg- 
ing rate evolution. It also reports for 30% of the current stellar 
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Mergers/ interacting phase: 
During the last 8 Gyrs, most lu- 
minous galaxies are expected to 
have experienced a major merger. 
During a major collision with a 
companion with mass above 1/4 
of the main galaxy mass, the disk 
is naturally suppressed as matter 
is falling to the mass barycenter, 
producing a compact morphology. 
Associated with short (0.1 Gyr) 
and strong peaks of star formation 
(most of them are LIRGs) 



Compact galaxy phase: 

Merger remnants linked to the for- 
mation of a significant fraction of 
stars in bulges : they show very 
blue central colours. Correspond 
to a decrease over 0.6-2 Gyr of 
the enhanced star formation due 
to the merging. Additional oc- 
currence of gas infall may subse- 
quently wrap around the bulge to 
form a new disk-like component as 
seen in many LCGs. 



Growth of disk phase: 

The star formation spreads over 
the entire forming disk as fed by 
large amounts of gas infall (from 
internal reservoir or accreted from 
the vicinity). Some of them show 
irregularities, which is interpreted 
to be related to minor merger 
events. Most distant spirals show 
a core much bluer indicating a 
recent and active star formation. 



Fig. 6. A sequence of z > 0.4 galaxies illustrating our proposed scenario for a recent formation of spirals in 3 major phases 
(major merger, compact galaxies, growing disks) 



Table 3. Characteristic times that spiral progenitors spend at each phase from z=l to z-0.4 {teiapsed - 3.3 Gyrs) 



Phase 


observed 


characteristic time 


characteristic time 


theoretical 






fraction 


(= fraction x 3.3 Gyr) 
(all spiral involved) 


(75% of spiral involved) 


expectations 


References 


Major merger (*) 


8% 


0.26 Gyr 


0.35 Gyr 


0.2-0.5 Gyr 


(1) 


Compact (major 


24.5% 


0.8 Gyr 


1.1 Gyr 


0.6-2 Gyrs 


(1), (2), (3) 


merger remnant) 
Irregulars (minor 


11.5% 


0.38 Gyr 


0.38 Gyrs 


0.2-0.5 Gyrs 


(1) 


mergers) 
Spirals 


56% 


1.8 Gyr 









- (*) corresponds to close pairs for which the two nuclei are separated by less than one galactic radius 

- (1): Tissera et al (2003); (2) Baugh et al (1996); (3) Cox et al (2004) 



mass formed over the last 8 Gyrs, which is essentially formed 
in gas rich spirals, those having approximately doubled their 
masses since z = 1. This is particularily consistent with the 
evolution of the L-Z relationship (Figure 4). The blue cores of 
distant spirals are also easily reproduced if most spirals were 
recently rebuilt through major mergers. The proposed evolu- 
tionary sequence is efficient in producing large bulges, and 
predicts a number density of late type spirals which sharply 



decreases from z - 1 to z = in agreement with the obser- 
vations (see Lilly et al. 1998). The relatively small fraction of 
present-day late type spirals (~ 20%) could then correspond 
to spirals which have not experienced a recent major merger. 
Although more qualitatively, the relative fraction of LCGs and 
major mergers selected from optical and IR wavelengths can 
be also interpreted: there are more optically selected LCGs 
than mergers (1st column of Table 2) because the merger rem- 
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nant time-scale is larger than that associated with the collision 
(see Table 3 and Tissera et al, 2002); however LIRGs include 
an almost equal number of LCGs as mergers (2nd column of 
Table 2) because of the lower star formation expected in the 
remnants. Finally a "spiral rebuilding" scenario is also consis- 
tent with an almost non-evolving spiral number density. Some 
major mergers -but not all- could involve 2 intermediate mass 
spirals, in which case a decrease of the spiral number density 
is expected. However this could be easily compensated for by 
merging of some smaller units whose end-products could enter 
the intermediate mass range. 

An important uncertainty is related to the actual merger rate 
since z=l, which requires accurate velocity measurements of 
both components. At z~ 0.3, Patton et al (2002) found that 
one fourth of interacting galaxies are unlikely to be merg- 
ing systems since their velocity difference exceeds 500 kms^^. 
Changes in the merger rate fraction would affect the relative 
fraction of (LCGs and LIRGs) events related to either major 
mergers or tidally triggered starbursts by galaxy interactions. 
This would not affect the crux of our scenario (Figure 6 still 
stands), i.e. an episodic star formation history for intermediate 
mass galaxies, which leads to major changes of their morpholo- 
gies as expected in hierarchical scenarios. If the major merger 
event per galaxy was significantly lower than 75%, it would 
imply that minor mergers and collisionless galaxy interactions 
are efficient enough to account for the strong and episodic star 
formation seen from IR observations. This "soft" alternative of 
our scenario requires a very (extremely?) large gas content in 
distant galaxies. 

Table 4 also tests a third scenario (called "collisional star- 
burst") assuming a marginal impact of major mergers. Such a 
scenario has been developped by Somerville, Primack & Faber 
(2001), and does not account for collisionless galaxy interac- 
tions. It succeeds in interpreting many galaxy evolution fea- 
tures up to z=3, but failed in reporting the star formation densi- 
ties derived from ISO and SCUBA experiments. It would also 
be difficult to reconcile it with the numerous populations of dis- 
tant compact galaxies and with the blue cores of most distant 
spirals. 

7. Discussion 

To firmly determine the above scenario will be very de- 
manding of observing time, requiring deep multi-wavelength 
observations (photometry & spectroscopy) for a much larger 
sample than the one presented here. Our results provide a good 
support for the hierarchical galaxy formation framework since 
most of the star formation is occurring episodically; at the 
level of our modest statistics, we predict that 42, 22 and 36% 
of the IR (episodic) star formation density is related to major 
mergers, minor mergers and gas infall, respectively. On the 
other hand, detailed simulations are also needed, especially to 
derive more accurate time-scales, and precise the mechanisms 
controlUng the disk growing (Birboim and Dekel, 2003; Dekel 
and Bimboim, 2004, in preparation). Given the fact that our 
scenario links distant and local galaxy properties in a simple 
way, we investigate here some of its potential failures. 



A major problem in reporting evolutionary trends of galaxy 
evolution from all inclusive surveys is related to Mahnquist 
biases, which severely limit their statistical significance. 
Indeed all inclusive surveys particularly are not adapted to 
follow galaxy evolution in the frame of hierarchical models, 
because the bright (generally massive) objects seen in the 
highest redshift bin are unlikely the progenitors of fainter (less 
massive) galaxies seen at lower redshifts. Here we simply 
consider all intermediate mass galaxies selected in a high-z 
volume (0.4 < z < 1) and then interpret their changes after 
relating them to local galaxies in the same mass range. Clearly 
our scenario prediction (most galaxies have recently doubled 
their masses) introduces a bias when comparing to local 
galaxies, but we find it produces a marginal effect, since the 
relative fractions of present-day E/SO and spirals are found 
barely dependent on changes in the low mass cut-off. 

One can also cast doubt on the fact that the CFRS I- 
selected sample is representative of all intermediate mass 
galaxies in a given redshift range. Our conclusions indeed 
account for a significant lack of red quiescent galaxies in the 
CFRS. Also, could we have missed a significant population 
of emission line galaxies ? Figure 7 shows the rest-frame 
colour-colour diagram for our sample. Independently of their 
nature (LIRGs, LCGs & normal starbursts), most galaxies 
are found within {{U - VOo = 1 ± 0.4 and (V - 1.4/im)o = 
1.3 + 0.5. There is only a hint that normal starbursts show 
bluer colours than LCGs and LIRGs, which is not unexpected 
if they are less affected by dust extinction. It means that a 
selection of galaxies using K band (rest-frame 1.4//m) or I 
band (rest-frame V) should find essentially the same mixture 
of emission line galaxies. Using the CFRS I-band limit we 
might have missed only a few dusty enshrouded starbursts 
which are similar to the LIRGs and LCGs at the top of Figure 7. 

Have we overestimated our stellar masses and misinter- 
preted a possible population of low mass galaxies in our sam- 
ple ? Stellar masses are estimated using various schemes for 
galaxy star formation histories, generally through a combina- 
tion of a (quiescent) underlying population with a burst (see 
Figure 7). Indeed our estimates are in the same range of those 
of BE2000 (same sample), agree within 30% with Franceschini 
et al. (2002) for LIRGs, and are in good agreement with those 
of Guzman et al. (2003) for LCGs (which correspond to their 
brightest sources in the HDF). However, the main uncertainty 
in the inferred stellar masses arises from the age of the stel- 
lar population, which leads us to suspect that masses of the 
most active star-forming systems (e.g LIRGs) possibly have 
been overestimated relatively to those of more quiescent galax- 
ies. This effect has been thoroughly investigated by Chariot 
(1998) who finds values of M/Lk as low as 0.2 for 0.1 Gyr 
old pure starbursts. LIRG spectral energy distributions (from 
UV to mid-IR and radio when available) are consistent with the 
superposition of a young and highly reddened burst to a stellar 
population found in local spirals (Flores et al. 1999). Assuming 
a 0.1 Gyr burst corresponding to the formation of 20% of the 
mass would lead to decrease the M/L^ ratio for LIRGs from 
~1 to 0.5, accounting also for the modest effect of extinction at 
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Table 4. Test of three scenarios 



Scenario/Evolution of 



"galactic downsizing" 
Intermediate mass 

galaxies already formed at z=l 
(BE2000, Cowie et al,1996) 



"spiral rebuilding" 
Minor/major mergers 
75% of Sp reprocessed 
(our scenario/this paper) 



"coUisional starbursts" 
Mostly minor merging 

(Somerville et al, 2001) 



References 



Mass growth 




+ 


+ 


H2004 










UlCKinson et (zUU4 ) 


L-Z (O/H) 


— 


+ 


+ 


Liang et al (2004b) 










this paper 


Pair statistics 


— 


+ 


— 


Le Fevre (2000) 










Bundy et al (2004) 


IR light evolution 




+ 




Flores et al (1999) 










Elbazetal(1999) 


# density of E/SO 


+ 


+ 


+ 


Schade et al (1999) 


# density of spirals 


+ 


+ 


+ 


Lilly et al (1998) 


# density of peculiar 


? 


+ 


? 


Brinchman (1998) 


(incl. LCGs) 








this paper 


spiral core colours 




+ 




Ellis et al (2001) 










Zheng et al (2004a,b) 



Note : a passing grade (+) means that the observables are consistent with the scenario, a negative grade (-) that there seems to be a significant 
inconsistency and a question mark indicates some doubt about the scenario consistency. 



that wavelength. The effect does not severely affect the fraction 
of LIRGs in intermediate mass galaxies (it drops from 15% to 
12%), and does not affect the overall mass production related 
to LIRGs. However the median doubling time-scale is then re- 
duced to 0.4 Gyr, leading to the formation of ~ 30% of the mass 
during a single infra-red event. A more rapid star formation in 
some galaxies with masses just below 3 lO'^M© would have 
also the advantage of more naturally relating distant galaxies 
to present-day intermediate mass galaxies. The detailed mod- 
elling of the stellar features of the VLT-FORS2 moderate reso- 
lution spectra of distant LIRGs will provide strong constraints 
on their past star formation histories, the intensity of the burst, 
its duration and the stellar mass fraction bom during the burst 
(Marcillac et al. 2004). 

We found that our results (morphological & spectral 
properties, LIRG fraction) are not significantly affected by 
adopting different schemes in assessing the stellar mass or 
by using different broad band selection criteria for our high-z 
sample. All Mb < -20 distant galaxies in the CFRS are found 
to be intermediate mass galaxies whether they are quiescent, 
"normal" starbursts, LCGs or LIRGs. In fact, we find that 
galaxy continuum properties (absorption line spectrum & 
broadband colours) are mostly unaffected by the burst strength, 
which can be simply understood if the burst brightening is 
somewhat balanced by the effect of dust extinction. In other 
words, optical properties of most LIRGs mimic those of normal 
starbursts. We conclude that selecting z< 1 galaxies from their 
apparent I to K band flux does not lead to a significant bias 
against a given class of galaxies, and that, to first order, what 
we see from these selections is what is there. This is at odds 
with the idea that present-day dwarf galaxies were playing 
a significant role in high z samples (Cowie et al. 1991; Koo 
et al. 1992) simply by a huge burst-brightening in optical. 



Besides the fact that this idea does not match with distant 
galaxy spectra and colours, another argument against the dwarf 
hypothesis is provided by the metal abundances of distant 
galaxies: although lower than typical values for intermediate 
mass spirals, these are much larger than expectations from 
progenitors of present-day dwarves. As an aside, one early 
motivation for our study of distant LCGs was to detect the 
[OIII]4363 Une which is ~l/3 the intensity of the Hy fine in 
local poor metal galaxies. No such detection (except for one 
source) has been made using deep exposure spectroscopy with 
the VLT, and the LCG spectra have conversely revealed a 
continuum with strong absorption lines, as expected from old 
stellar component. 

An important assumption of our scenario is the requested 
availability of supplementary gas to progressively form the 
stellar disk content after the compact phase. Much higher gas 
infall was expected in the past as shown by the large redshift 
increase of the number densities of starbursts and LIRGs. A 
high gas reservoir in progenitors of spirals is not unexpected 
at these epochs, because 10 Gyr ago (z = 2), most of the 
baryonic mass was in the form of ionised gas, whose mass 
content has by now dramatically decreased, by a factor ~ 50, 
at z = (Fukugita, Hogan and Peebles, 1998), even when 
accounting for neutral and molecular gas. We argue that this 
change of state of the baryonic matter did not take place only 
between z = 2 and z =1 (elapsed time 2.6 Gyr), but is likely 
the result of a longer process over the period between z = 2 
and z = 0.4 (elapsed time 6 Gyr). Cosmic chemical evolution 
models (Pei, Fall and Hauser, 1999) successfully reproduce 
the CSFD evolution, and assume a decrease of the ionised gas 
density by a factor 5 since z = 1. Indeed in the distant large 
spirals that are LIRGs, the star formation process is spread aU 
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1 2 

(U - V)o (AB) 



Fig. 7. Rest-frame color (V - 1 Alb/umjAB versus (U - V)ab for CFRS luminous galaxies (Mb < -20, same symbols as in Figure 
4). It is adapted from Hammer et al. (2001, see their Figure 5), who have shown that the rest-frame colours used here match 
particularily well the observed colours in the redshift range z=0.4-0.8, leading to very small errors related to k-corrections. Here 
we represent CFRS galaxies for which the photometric error bars on colours are below 0.4 magnitude. Lines show Bruzual and 
Chariot (1999) models with solar metallicity and an exponential declining SFR with an e-folding time-scale t. The dotted line 
show the model with t = 0.1 Gyr, the short-dashed line the model with r = 1 Gyr and the long-dashed line the model with a 
constant SFR. The bold short-dashed line represents the t= 1 model with an extinction Ay - I applied. The bold dot-short dashed 
line represents the t - I Gyr model, assuming a metallicity Z/Zq -3.2. A composite model, where 80% of the light is contributed 
by an old stellar population formed 5 Gyr ago in an instantaneous burst, and the remaining 20% of the light by a younger stellar 
population (from t=0 to t=lGyr) with t-1 Gyr, is shown as a solid (blue, almost horizontal) line. Red (almost vertical) line show 
the superposition of a normal starburst color (cross) with an extincted starburst (see text). A small arrow on the upper-right side 
of the Bruzual and Chariot (1999) models indicates the uncertainties in the modeling of the old stellar population, as discussed 
by Chariot et al. (1996). 



over the disk, as shown by its blue colour (Zheng et al. 2004a), 
as expected if large amounts of gas are being accreted. 



8. Conclusion 

We have re-examined here the properties of distant galaxies 
mostly on the basis of a panchromatic follow-up study of the 
CFRS. We find that: 

- SFRs based on the [01I]/13727 emission line lead to 
severely underestimated values, which can lead to mislead- 
ing results on the star formation history of distant galaxies 
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- 15% of intermediate mass galaxies at 0.4< z < 1 are in- 
deed LIRGs and the associated star formation density is 
sufficient in itself to form 38% of the total stellar mass in 
present-day intermediate mass galaxies 

- all emission Une galaxies (starbursts and LIRGs) at z~ 0.7 
have gas phases which are metal deficient by a factor 2 
when compared to those of local spirals 

- most distant field galaxies (but E/SO) show significant star 
formation in their Ikpc radius central cores 

- optical/spectral properties of LIRGs are so similar to those 
of other galaxies that only IR measurements are able to de- 
scribe how the star formation density is distributed between 
the different galaxy types 

- the only way to interpret the high occurrence of LIRGs is 
to associate them to episodic star formation histories for 
most galaxies, in agreement with the hierarchical galaxy 
formation scenario 

These episodic and violent star formation events can be 
simply associated to major changes in galaxy morphologies. 
We have then tested three galaxy formation scenarios against 
eight rather robust features of galaxy evolution. The best 
scenario accounts for a recent formation of bulges in spirals 
through merging or tidally triggered central bursts, followed 
by an inside out formation of the stars in disks. In this scenario, 
at least half of the bulge stellar content was made earlier in 
their progenitors, before the last major phase of accretion. It 
then provides a simple origin for the reported recent stellar 
mass production (30 to 50% of the stellar mass formed since 
z=l) and for a recent peak in the cosmic star formation history 
(H2004). It is consistent with the simultaneous decreases, 
during the last 8 Gyrs, of the UV and IR luminosity density, 
of the merging rate, and of the number densities of LIRGs and 
of LCGs. It also naturally accounts for all the major changes 
of the galaxy morphologies in the last 8 Gyrs, as well as for 
an approximately unchanged density of large disks at z ~ 1. 
Because of the striking similarities between properties of 
galaxies at different luminosities (average 0/H abundances, 
fraction of merging events and of compact galaxies), possibly 
most present-day spirals, independently of their mass range, 
have recently experienced their last major phase of the building 
of both their bulges and their disks. More than a third of the 
present-day stellar mass formed at z < 1 is associated with 
more than half of the stellar mass formed in spirals during 
violent events which have shaped them as the Hubble sequence 
describes them today. 

All the important numbers used here to estabhsh our 
scenario (e.g. fraction of LIRGs and of various morphological 
types, SFR, stellar masses etc..) have been systematically 
cross-checked and found in good agreement with the littera- 
ture. We also find that biases related to the I band selection 
criterion used in the CFRS provide minor effects. This is 
likely due to the fact that dust-enshrouded starbursts (such 
as LIRGs) show optical properties which mimic those of 
normal starbursts, in such a way that optical observations are 
generally unable to distinguish them. This is well reported 
in our scenario: the addition of a highly reddened burst to 



a normal high-z galaxy has no severe impact on its optical 
continuum. The way IR Ught impacts our optical view of the 
galaxy evolution is then rather subtle and it explains well the 
often associated wording "hidden star formation". Without 
observing in the IR and accounting only from the UV Ught, 
all distant galaxies apparently show modest star formation rate 
and stellar mass increase. 

A significant number of recent major mergers is qualita- 
tively consistent with the dramatic increase of type 1 and 2 
AGN out to z=l reported by Miyaji, Hasinger & Schmidt 
(2000). Major mergers are particularily efficient in forming 
large bulges which can account for the numerous population 
of present-day early-type spirals (see also Zheng et al, 2004b). 
The link between distant and nearby galaxies could be fur- 
ther investigated: what is the relative impact of both star for- 
mation and mass aggregation at z< 1 for present-day spirals 
and ellipticals ? Could we match our scenario with the Milky 
Way and M3 1 past histories ? While it is unlikely that a ma- 
jor merger has recently occurred in the Milky Way, the M31 
system appears to be much more complex and the significant 
intermediate age population in its halo supports the idea of a 
recent galaxy merger (Brown et al, 2003). It is also necessary to 
predict more accurately the relative fraction of star formation 
occurring during major mergers, minor mergers, interactions 
and gas infall, respectively. This requires unambiguous esti- 
mates of the actual merging rate and of the fraction of LCGs 
which are end products of merging events. We believe that 
estimates of velocity fields of distant galaxies with the VLT 
(with FLAMES/GIRAFFE) will address these issues, as well 
as the evolution of the Tully Fisher relationship. Large follow 
up of the deepest fields studied by HST and Spitzer should be 
done, using moderate spectral resolution spectroscopy to derive 
SFR, extinction and chemical gas properties, as well as velocity 
fields provided by 3D spectroscopy. 
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